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Abstract 

An intrinsic gap depression at the Superconductor-Insulator interface due 
to the very short value of the coherence length in High-Tc Superconductors 
[HTSs] is considered, in the framework of a mixed (s+id)-waye pair symme- 
try for the order parameter ranging from pure s to pure d-wave. This gap 
depression acts as the main physical agent causing the relevant reduction of 
Ic (T) Rn{T) values with respect to BCS expectations in HTS SIS Josephson 
junctions. Good agreement with various experimental data is obtained with 
both pure s-wave and pure d-wave symmetries of the order parameter, but 
with amounts of gap depression depending on the pair symmetry adopted. 
Regardless of the pair symmetry considered, these results prove the impor- 
tance of the surface order-parameter depression in the correct interpretation 
of the Ic{T)Rn{T) data in HTS SIS junctions. In a case of planar YBCO- 
based junction the use of the de Gennes condition allowed us to tentatively 

obtain an upper limit for the amount of d-wave present in the gap of YBCO. 
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As a consequence of the very short coherence length ^ of high-Tc superconductors [HTSs], 
an intrinsic depression of the order-parameter at Superconductor-Insulator [S-I] interfaces 
arises which can represent one of the main reason for the reduced Ic{T)Rn{T) values in 
HTS SIS 0,1 and SIS' Josephson junctions §. 

Also pair symmetry is expected to contribute to such experimental results, and, in the 
present paper, a model is developed by us which takes into consideration pure s-wave, pure 
(i-wave and mixed (s + 2(i)-wave pair symmetry together with a suitable gap depression, 
extending and generalizing our previous model developed for T > only in the pure s-wave 
case 

By making use of the space-dependent expression of the gap A(x, T)= A (T) 
tanh [(x + xo)/ (V2C(T))] • 9 {x — w) which may be derived from Ginzburg-Landau equa- 
tions, Ic{T)Rn(T) values in SIS and SIS' Josephson junctions [JJs] may be evaluated. Here 
2w is the thickness of the insulating layer and xq accounts for the spatial slope of A (x, T) 
near interfaces 0. As in Ref. we assume that the d-wave component of the order pa- 
rameter Arf has the same spatial behaviour as the s-wave component A^ and, moreover, 
we neglect the fourth-order terms in the Ginzburg-Landau equations for the (s + id)-wave 
pair symmetry. We imagine to " section" the gap A (x, T) into independent channels SAi 
IP,^, each of them giving rise to a parallel contribution 61^ (T) and ^G^. (T) to the total 
critical current and normal conductance, respectively. The physical situation just described 
is depicted in Fig. 1 where our model for a HTS SIS Josephson junction is shown. By 
summing and averaging over all the parallel contributions, that in general differ since each 
"slice" corresponds to a different thickness of the barrier (see Fig. 1), we obtain: 

/c(T)^^E[4 (T)-5A, (x,T)] 
Gn{T)c^^^Y.[Gn. {T).SA, (x,r)] 



In the continuous limit these expressions become: 
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and 



[Gn (x, T)]-^ = Rn (x, T) = i?o (T) ■ 10'(^-^) 

where Rq is a parameter representing the resistance of the channel at x = w, which 
gets ehminated after performing the integral for Gn (T), and where the space dependence 
of Rn (x,T) has been heuristically derived under global conditions of consistency The 
expression for the spatial dependence of Ic (T) in mixed (s-/-zd)-wave symmetry is derived 
by introducing the spatial dependencies of A,, and A^ in the results of Ref. : 
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where ui = {21 — ViT^ksT is the Matsubara frequency and K is the complete elliptic 
integral of the first kind. 

After performing rather cumbersome calculations, the general expression for the temper- 
ature dependent critical voltage /c(T)_R„(T) in the case of SIS JJs, in mixed (s-/-zc?)-wave 
symmetry and in presence of surface depression of the order parameter is given by: 
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KJl,T,z) = K 
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and . (X, T) = tanh (^) , 7 (T) = tanh (^) , ^ (T) = ^ (T) (i^) . 

In the previous expressions we considered A^ (T) = [1 — e) ■ Agx{T) and Ad{T) = 
e Aex (T) where e is the fraction of c?-wave present in the order parameter, Agx is the 
experimental value of the gap determined, for example, in tunneling experiments []^,^ and 
where we have used for the temperature dependence of A (T) and ^ (T) the standard BCS 
expressions. In the special case of pure s-wave pair symmetry Eq. (2) reduces to 0: 

k(i^)%anh(^5^) + Li;^(l-^)%anh(^5^) dz 

[J-cK-L )rin\-L )\sis — ^ /i \i9 1 /, \i? 



which tends to the Ambegaokar-Baratoff [AB] model |^ in the limits valid for a low-Tc 
superconductor (LTS) (i.e. xq — > +00,,^ (0) +00, xq/^ (0) +00): 

lim \IAT)RJT)]Z^ s+dep.gav ^ ^A, (T) / A^\ ^ 

HTS^LTS^ ^ ' ^ ^J^^'^ 2e \2kBT J 

When gap depression at S-I interfaces is neglected, the result obtained by Xu et al. 
in the general mixed {s + Z(i)-wave case is reproduced: 
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On the other hand, in the case of pure c?-wave pair symmetry and in presence of a 
depression of the order parameter all the components of Eq. (2) reduce to the analogous 
expressions for A^iT) = (i.e. e = 1) : 

[/,(T)i?„(T)]^7;'+'''p-^'^^ = lim [/,(T)i^„(T)](f+;')+''^^•^"^ 



We remark that the free parameters of the model expressed by Eq. (2) in the most 
general case are three: w, Xq and e. When the thickness of the insulating barrier is known 
(for example in planar SIS JJs) the number of free parameters reduces to two. It becomes 
only one if, in a planar SIS junction, the amount of d-wave is also fixed as in the pure s-wave 
case {e = 0) or in the pure d-wave one (e = 1). 

A reasonable estimation of the minimum (intrinsic) amount of gap depression present at 
S-I interfaces may be obtained by calculating the xq value - from now on called Xq^ - that 
derives from the de Gennes condition at the S-I interface [|T0[] : 
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where, in the hypothesis of Ginzburg-Landau behaviour for A (x), we have 

1 _ V2 1 

The previous expression, together with the approximate value & — (0) /{a — w) deter- 
mined by de Gennes, where a is the range of the x axis where A (x) varies appreciably [jlO 



permits to calculate xf^. When, due to a short (0), a large intrinsic gap depression at S-I 
interfaces is expected, ^^^^ < 0.4 so that b ~ w + Xq'^. This condition is always satisfied in 
the HTS junctions we will describe thereafter. If we impose that A (a, T) /A (oo, T) = 0.99, 
considering that 0.99 ~ tanh2.6, a value of a may be yielded so that an estimate for Xq*^ is 
obtained by means of Xq'^ ~ 0.6 • ,^ (0) — w. In principle, the xq values determined by fitting 
the experimental data with our depressed-gap model have to be comprised between Xg*^ and 
Xq = —w that corresponds to a complete gap depression at the interface i.e. A {w)= 0. 

Agreement of such a model with experimental Ic{T)Rn{T) data in HTS SIS JJs has 
been tested for the two cases of pure s-wave and pure d-wave pair symmetry. On the other 
hand, a fit of the experimental curves by the complete (s-/-zc?)-wave model in absence of 
indications concerning the amplitude of the d-wave component of the gap would have been 
of very modest interest due to the numerical impossibility to separate the contributions of 
the gap depression xq and of the d-wave gap e to the lowering of Ic(T)Rn(T). As we will 
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see in the following examples, results are very good with both types of pair symmetry, but 
with a different amount of gap depression in each case. 

Examples presented in this work include one of our recent Ic{T)Rn{T) behaviours 
obtained from reproducible nonhysteretic I-V curves in Bi2Sr2CaCu2 08+a; (Tc = 85 K, 
Abscco (0) = 23 meV, ^bscco (0) = 1-67 nm) single-crystal high-quality Josephson break 
junctions (Fig. 2). RSJ model has been used to obtain Ic{T) and Rn(T) values as these 



latter were high enough to permit to treat the junctions as tunnel SIS ones 13,0] . Best-fit 
values for the amount of the gap depression are (xq)^ = —1.75 nm and (xq)^ = 0.5 nm, 
while the thickness of the barrier, fitted by using the pure s-wave expression and held con- 
stant also for the pure d-wave case, is 2w = 4.4 nm. Figure 2 shows the results of the 
fit of experimental data by using the depressed-gap model compared to the standard BCS 
results (no gap depression) in s- and d-wave symmetry. As in all the figures from 2 to 5, 
open symbols are experimental Ic(T)Rn(T) data, dash-dots and dots represent predictions 
obtained considering pure s-wave (AB model 0) and pure d-wave pair symmetry without 
any gap depression, respectively, while the results of our model applied to the pure s-wave 
and the pure d-wave cases are reported with a dashed and continuous line, respectively. In 
all the presented cases the superconducting gap A(0) has been determined from tunneling 
experiments or from literature while the coherence length C, (0) was taken from literature. 

The depressed-gap model has been tentatively used also in Bi2Sr2Cu06+x (Bi-2201) 
single-crystal break-junctions [|13| where Tc = 9.6 K, Absco (0) = 3 meV and ^bsco (0) = 4 



nm. The results of the fits in these nonhysteretic Bi-2201 break-junctions are shown in Fig. 
3 and the fitted values of the gap depression are equal to (xq)^ = —4.0 nm, (xq)^ = 0.5 nm. 
The barrier thickness, as usual determined by the s-wave depressed-gap fit, is 2w = 9.8 nm. 
Fig. 4 reports data taken from a SIS YBCO-based bicrystal grain-boundary junction 



fabricated on a SrTiOs (STO) substrate |T^ and fitted by our model with a barrier thickness, 
determined as in the previous cases, 2w = 2.0 nm . Here the amounts of gap depression 
are (xq)^ = —0.6 nm and (xq)^ = 1.0 nm, while = 85 K, Aybco (0) = 16 meV and 
^YBCO (0) = 2.8 nm. 
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Fig. 5 reports the case of a YBCO/STO/YBCO multilevel edge J J |T^J^ with ex situ 



oxygen plasma treatment, which is well fitted by our model (both s- and c?-wave symmetry) 
with only one adjustable parameter, being 2w = 7.2 nm already known independently from 
the fabrication process. The fitted values of the parameter xq are (xo)^ = —3.3 nm and 
(xq)^ = —0.4 nm. The physical values characteristic of the superconducting state used in 
the fit are: Tc = 85 K, Aybco (0) = 16 meV and ^ybco (0) = 2.8 nm. It is clear from Figs. 
2 to 5 that both the s-wave and the d-wave depressed-gap curves fit properly the magnitude 
and, in most of the cases, also the shape of the various experimental Ic{T)Rn{T) data, 
but with different amounts of gap depression: pure (i-wave fit requires a smaller amount of 
surface gap depression (i.e. a larger value of Xq) than pure s-wave one. 

In the latter junction, taking advantage of the a priori knowledge of the exact value of 
w, the behaviour of the amount of gap depression (accounted for by Xq) vs. e (percentage 
of d wave) has been obtained by using the general mixed (s + i(i)-wave model of Eq. (2) 
and reported in Fig. 6. Since in that case x^^ ^ —1.9 nm, from the graph of Fig. 6 a pure 
(i-wave would be excluded, and an upper limit of the order of 40% for the percentage of 
(i-wave pair symmetry in YBCO is draftily given. Anyway we must remark that the results 
are very sensitive to the gap value and, therefore, different values for Aybco (such as 20 
meV rather than the 16 meV we used following the authors of Ref. ||15|) would yield for xq 
values compatible also with a pure d-wave pair symmetry. The Xq*^ values, calculated by 
means of the above mentioned approximate expression for the junctions of Fig. 2, 3 and 4, 
are always comprised between (xo)<i and (xo)d- Even considering the mentioned sensitivity 
of the results to the value of the gap and to the value of w, this fact seems to indirectly 
question the presence of pure d-wave symmetry in the HTS materials of these junctions. 

With this extension of our model, any comparison with the experimental results of SIS' 
JJs |T6|JT7[1 (where S' is a low-Tc superconductor) is no longer meaningful due to the orthog- 
onality between wave functions with symmetries s and d. Actually, the mere detection of a 
Josephson current in SIS' JJs rules out the existence of pure d-wave pair symmetry. Only 
depressed-gap pure s-wave symmetry |Q or mixed (s-/-z(i)-wave pair symmetry plus order- 



parameter depression can be acceptable in SIS' JJs. In the latter case, if the ratio A^/Ag 
is not fixed (i.e., if it's a free parameter of the model), its effect in reducing /c(T)i?„(T) 
values becomes not distinguishable and not separable from the effect eventually operated by 
a S-I gap depression. On the other hand, considering the ratio A^/A^ as a free, adjustable 
parameter leads to an unphysical sample-dependence of the pair symmetry which does not 
appear acceptable. 

A first conclusion is therefore the full confirmation of the important role played by a 
surface gap depression in reducing the "quality factor" Ic{T)Rn{T) below its theoretical BCS 
value (regardless of the type of gap symmetry considered) in HTS SIS Josephson junctions. 
As it is shown in Figs. 2 to 5 and in the framework of our model there is no possibility 
to fit the shape and, especially, the magnitude of experimental data by using pure 5-wave, 
pure d-wave or mixed {s+id)-wsive models without any gap depression. Moreover, we remark 
that fitting of reduced Ic (T) /Ic (0) values rather than Ic (T) ones is not meaningful - even 
if very often present in literature - since the main discrepancy of experimental results from 
theoretical predictions is in the magnitude rather than in the shape of the curves. 

As we have shown in the previous examples, the value of Xq accounting for the depression 
of the order parameter depends on the type of pair symmetry adopted and, therefore, we can 
also conclude that - when applied to several experimental results - this model might yield a 
contribution to the pair symmetry debate only provided the amount of gap depression at S-I 
interface in each case examined is given independently from other measurements or physical 
considerations, and is not a free, adjustable parameter as we have treated it in our present 
calculations. 
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FIGURES 



FIG. 1. Surface-depressed order parameter A (x, T) as function of x and a "slice" SA^ 
into which we imagine to section the gap, in order to compute Ic(T)Rn(T) in SIS JJs. 

FIG. 2. Comparison of the experimental Ic{T)Rn{T) obtained from a Bi2Sr2CaCu208+a; 
single-crystal Josephson break junction with the curves yielded by various theoretical models 
(see text): open symbols are experimental data; dash-dots and dots represent predictions 
obtained considering pure s-wave and pure d-wave pair symmetry without gap depression, 
respectively, while the results of our model both in the pure s-wave and in the pure d-wave 
cases are reported with a dashed and continuous line, respectively. 

FIG. 3. Comparison of the Ic(T)Rn{T) behaviour obtained from a Bi2Sr2Cu06+x single- 
crystal break-junction with the results of various models as shown in Fig. 2. 

FIG. 4. The same as in Figs. 2 and 3 but with the /c(T)i?„(T) data taken from a 
YBCO-based bicrystal grain-boundary junction fabricated on SrTiOs substrates (from Ref. 

0). 

FIG. 5. The same as in Figs. 2-4 but with the experimental /c(T)i?„(T) data taken 
from a SIS YBCO/STO/YBCO multilevel edge JJ (from Ref. [|T§) and fitted by our model 
with only one adjustable parameter. 

FIG. 6. Behaviour of the amount of gap depression (accounted for by Xq) vs. percentage 
of d wave (in the general mixed (s -|- i(i)-wave model), estimated taking advantage of the a 
priori knowledge of the barrier thickness 2w in the case of Fig. 5. 
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